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We report the observation of the X(3872) in the J/ipir + n~ channel, with J/tp decaying to n + fj,~ , 
in pp collisions at y/s= 1.96 TeV. Using approximately 230 pb~* of data collected with the Run II 
D0 detector, we observe 522 ± 100 X(3872) candidates. The mass difference between the X(3872) 
state and the J/tp is measured to be 774.9 ± 3.1 (stat) ± 3.0 (syst) MeV/c 2 . We have investigated 
the production and decay characteristics of the X(3872), and find them to be similar to those of 
the ip(2S) state. 

PACS numbers: 14.40.Gx, 14.40.Lb, 13.25.Gv, 13.25.Jx, 13.85.Ni 



A new particle, the X(3872), was recently discov- 
ered by the Belle Collaboration in the decay mode 
B ± -» XK ± (X -» J/tpn+n-) where the mass of the 
X(3872) was measured to be 3872.0 ± 0.6 (stat) ± 0.5 
(syst) MeV/c 2 . The existence of the A(3872) state has 
been confirmed (also decaying in the J/ipir + ir~ mode) in 



pp collisions by the CDF Collaboration Q. At this time, 
it is still unclear whether this particle is a cc state, or a 
more complex object. See for example 

The charmonium state ip(2S), with mass rn^(o,s) = 
3685.96 ± 0.09 MeV/c 2 0, has the same decay mode, 
providing a good benchmark for comparison with the 
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A(3872). The ipi^S) mesons produced in pp collisions 
can originate either from decays of B hadrons or from 
direct production. The tp(2S) mesons from B decays 
have longer effective decay lengths and tend to be less 
isolated than directly produced tp(2S) mesons 0- 

We examine the production rate of the A (3872) rela- 
tive to ip{2S) as a function of the transverse momentum 
with respect to the beam axis (pt), isolation and decay 
length, as well as a function of rapidity {y — \ log j^p^ , 
where E is the energy and Pl is the longitudinal mo- 
mentum with respect to the beam axis), to determine 
whether the production characteristics of the X(3872) 
are similar to those of the ip(2S). We also compare the 
angular decay distributions of the 7r + 7r~ and sys- 
tems in X(3872) decays with those from ip(2S), to check 
for any differences in helicities of these two states. 

The data set used in this Letter was collected in pp 
collisions at ^=1.96 TeV between April 2002 and Jan- 
uary 2004, and corresponds to an integrated luminosity of 
approximately 230 pb -1 . The D0 detector is described 
elsewhere 9J. The components most important to this 
analysis include the vertex, central tracking and muon 
systems. The D0 tracking system consists of a silicon mi- 
crostrip tracker (SMT) and a central fiber tracker (CFT), 
both within a 2 T solenoidal magnetic field. 

The SMT has approximately 800, 000 individual strips, 
with typical pitch of 50-80 /an, and a design optimized 
for tracking and vertexing over the range \rj\ < 3, where 
f] = — ln[tan(0/2)] is the pseudorapidity and 6 is the po- 
lar angle measured relative to the proton beam direction. 
The system has a six-barrel longitudinal structure, each 
with a set of four layers arranged axially around the beam 
pipe, and interspersed with 16 radial disks. The system 
provides a resolution, in the plane transverse to the beam 
axis, for the distance of closest approach of a charged par- 
ticle relative to the primary vertex of ss 50 /xm for tracks 
with pt ~ I GeV/c, improving asymptotically to 15 /im 
for tracks with px > 10 GeV/c. 

The CFT comprises eight thin coaxial barrels, each 
supporting two doublets of overlapping scintillating fibers 
of 0.835 mm diameter, one doublet being parallel to the 
collision axis, and the other alternating by ±3° to provide 
information along the beam axis. 

The muon system is located outside the calorimeters, 
and consists of a layer of tracking detectors and scintilla- 
tion trigger counters in front of 1.8 T toroidal magnets, 
followed by two similar layers behind the toroids. Track- 
ing in the muon system in the range \rj\ < 1 relies on 
10 cm wide drift tubes |10|. while 1 cm mini-drift tubes 
are used for 1 < \rj\ < 2. 

J/ijj — > fj, + [i~ decays are selected by triggering on 
dimuons using a three-tier trigger system. The first trig- 
ger level uses hardware to form roads defined by hits 
in two layers of the muon scintillator system. The sec- 
ond trigger level uses digital signal processors to form 
track stubs defined by hits in the muon drift-chamber 
and muon scintillator systems. The third level comprises 
a farm of computer processors with access to the entire 



event. Events passing the third-level trigger are recorded 
for analysis. 
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FIG. 1: AM = M(n + fi~TT + TT')-M(p + fi') for all candidates 
satisfying the selection requirements. The solid curve is a fit 
to the data and the dashed curve represents the background 
under each peak. The insert shows the mass distribution of 
the J/tp candidates used in the analysis. 

Muons are identified by extrapolating charged particle 
tracks from the central tracking system that match with 
muon track segments formed from hits in the muon sys- 
tem. Oppositely charged muons are combined to form 
J 'ftp candidates, which are then combined with two op- 
positely charged particles assumed to be pions. At least 
two of these four tracks are required to have at least 
one hit in the SMT. To reduce background from com- 
binatorics, events are required to satisfy the following 
selection criteria. An event is required to have less than 
100 tracks. Jji\) candidates are selected by requiring the 
invariant mass of the fi + fi~ system to be between 2.80 
and 3.35 GeV/c 2 , and the transverse momentum with 
respect to the beam axis (pt) of the J/ip is required to 
be greater than 5 GeV/c. In addition, the px of each of 
the two pions must be greater than 0.7 GeV/c, and the 
spatial separation, AR, between the momentum vector 
of the J/ipir + ir~ system and each pion momentum vec- 
tor is required to have AR < 0.4, where AR is defined as 
y/ {A(f>) 2 + (Arj) 2 , with <j> being the azimuthal angle. The 
invariant mass of the two pions, M(7r + 7r~), is required 
to be greater than 0.52 GeV/c 2 , and the \ 2 01 a nt to 



the ^ + /x 



7T 7T 



vertex is required to be less than 16 (for 



five degrees of freedom). 

Figure ^ shows the distribution in the mass difference 
AM = A/(^ + ^i~7r + 7r~) - M([i + fi~), after all selections. 
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Superimposed is a fit to the data, where Gaussians are 
used to describe the ip(2S) and the A(3872), and a third- 
order polynomial is used to account for the background. 
The fitted width of the X(3872) peak is 17 ± 3 MeV/c 2 , 
which is consistent with detector resolution. The re- 
sults of the Gaussian fit yield 522 ± 100 A (3872) can- 
didates with AM = 771.9 ± 3.1 (stat) MeV/c 2 . The 
yield of the X(3872) is dependent on the value of the 
M(7r + 7r _ ) requirement. When we tighten the require- 
ment on M(ir + ir~) to be greater than 0.60 GeV/c 2 , 
0.65 GeV/c 2 and 0.70 GeV/c 2 and fix the mean of the 
A(3872) to its nominal value, we obtain yields of 433±98, 
336±78 and 170±54 events, respectively. 

The position of the J/tp mass in Fig. \l\is shifted rela- 
tive to the currently accepted value Q • This is within the 
uncertainties of our momentum scale, but we correct the 
measured AM by the ratio of the currently accepted mass 
difference between the ip{2S) and the J/ip @ and our cor- 
responding measured mass difference. We assign a 100% 
systematic uncertainty to this mass scale correction, giv- 
ing a final measurement of AM = 774.9 ± 3.1 (stat) ± 
3.0 (syst) MeV/c 2 . We investigated systematic effects 
by using different background parameterizations and bin 
widths in our fits, and found the largest change in the 
central value of the mass difference to be 0.6 MeV/c 2 . 

To investigate the characteristics of the X(3872) state, 
we study its production and decay properties, and com- 
pare the signal yields of the A(3872) to the ip(2S). For 
example, Fig. [21 shows the AM distribution in two sepa- 
rate regions of rapidity of the A (3872). The data is also 
separated into regions of effective proper decay length, 
transverse momentum of the X(3872) candidate, isola- 
tion, and decay angle. 

The effective proper decay length, dl, is defined as the 
distance in the transverse plane from the primary vertex 
to the decay vertex of the J / ip scaled by the mass of the 

system 
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7r + 7r system divided by the fi + fJ. i 
Pt- 

The isolation is defined as the ratio of the X(3872) mo- 
mentum to the sum of the momentum of the X(3872) and 
the momenta of all other reconstructed charged particles 
within a cone of radius AR = 0.5, about the direction of 
the X(3872) momentum. 

The helicity of the 7r + 7r~ (/j + /j~) system can be in- 
ferred by boosting one of the pions (or muons) and the 
X(3872) into the dipion (dimuon) rest frame, and mea- 
suring the angle 9^(9^) between them. The cosine of this 
angle is used for the comparison between the X(3872) 
and ?/)(2S). 

The widths of the ip(2S) and X(3872) are fixed to the 
values obtained in the fit to the full sample, and the 
fitted number of ip(2S) and X(3872) candidates in each 
region of chosen variables is given in Table [I] The data 
in all regions are well represented by our choice of fitting 
function, with fit probabilities all greater than 15%. In 
Fig. |3| we present the results of Table U in a graphical 





0.6 0.7 0.8 0.9 A 

"W - M^. (GeV/c 2 ) 



FIG. 2: AM distribution for two ranges of rapidity of the 
fi + 7r + 7r~ system. 



form. For example, we compare the event yield fraction: 
X(X(3872)) \y\ < 1 



/x(3872) 



to that for ip{2S) 



/■0(2S) 



N(X(3872)) |y| < 2' 



N(j>(2S)) jyj < 1 
N(j>(2S)) \y\<2- 



From the ratios of isolation and decay lengths, it is ap- 
parent that the production of the X(3872) has a similar 
mixture of prompt and long-lived contributions as the 
V(2S). 

In summary, we observe the X(3872) particle in ap- 
proximately 230 pb _1 of data with a significance of 5.2 
standard deviations. The mass difference between the 
X(3872) particle and the J/tp is AM = 774.9 ±3.1 (stat) 
± 3.0 (syst) MeV/c 2 . When the data are separated ac- 
cording to production and decay variables, we find no 
significant differences between the X(3872) and the cc 
state ip(2S). 
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TABLE I: Number of tp(2S) and X(3872) candidates for 
different ranges of variables. The fitted widths of the X(3872) 
and ip(2S) are not constrained in the initial selection, but in 
the different ranges they are fixed to the values obtained from 
the full sample. Apart from the full sample, uncertainties are 
only from the normalization of the fitted Gaussian function 
and do not include small contributions from uncertainties in 
the background. 



Reg 


;ions 




Number 
of V(2S) 


Number 
of A (3872) 


Initial selection 




1192 ± 55 


522 ± 100 


(a) 




) > 15 GeV/c 


396 ± 26 


179 ± 39 




pt(j yy> tt + ty~ 


) < 15 GeV/c 


796 ± 41 


358 ± 64 


(b) 


\y\ < 1 




741 ± 37 


316 ± 57 




1< M <2 




449 ± 31 


204 ± 49 


(c) 


cos(^) < 0.4 




589 ± 34 


288 ± 53 




cos(6^) > 0.4 




606 ± 34 


244 ± 53 


(d) 


dl < 0.01 cm 




838 ± 41 


351 ± 66 




dl > 0.01 cm 




359 ± 26 


164 ± 41 


(e) 


isolation = 1 




257 ± 20 


85 ± 22 




isolation < 1 




942 ± 44 


438 ± 72 


(f) 


cos(6» M ) < 0.4 




593 ± 33 


232 ± 46 




cos^) > 0.4 




602 ± 35 


288 ± 60 
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FIG. 3: Comparison of event-yield fractions for A(3872) and 
i)(2S) in the regions: (a) p T (J/ipTv + Tv~) > 15 GeV/c; (b) \y\ 
of J/ip iy + tt~ < 1; (c) cos(S^) < 0.4; (d) effective proper decay 
length, dl < 0.01 cm; (e) isolation = 1; (f) cos(6< M ) < 0.4. 
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